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Abstract 



We discuss calculations of three-body observables for the breakup of 8 B on 
a 58 Ni target at low energy using the coupled discretised continuum chan- 
nels approach. Calculations of both the angular distribution of the 7 Be frag- 
ments and their energy distributions are compared with those measured at 
several laboratory angles. In these observables there is interference between 
the breakup amplitudes from different spin-parity excitations of the projec- 
tile. The resulting angle and the energy distributions reveal the importance 
of the higher-order continuum state couplings for an understanding of the 
measurements. 
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I. INTRODUCTION 



Projectile breakup is an important reaction channel in the scattering of weakly bound 
nuclei. An accurate treatment of breakup is therefore a major ingredient in attempts to un- 
derstand the properties of light exotic radioactive nuclei from reaction studies. The number 
of published experimental breakup studies, and also their accuracy, has increased rapidly. 
These include reactions in which both nuclear and Coulomb breakup effects are expected to 
be significant, e.g. p1-|T7j|. Until recently, the low intensities of available rare isotope beams 
has meant that many of the experiments were either designed to measure inclusive cross 
sections, with incomplete kinematics, or did not have adequate statistics to allow the ex- 
traction of exclusive observables. The cross sections extracted from the measurements were 
often integrated over fragment energies or angles, or both, and inevitably some details of 
the physical process were lost as a result. This is no longer the situation. Secondary beam 
intensities are becoming sufficiently high that coincidence experiments are now practical 
and, in the future, data will more routinely require a fully three- or more-body study, e.g. 
T8| , p!9|| . The need for precise theoretical predictions of the breakup of two-body projectiles, 



and of their three-body observables, is the primary motivation for this work. 

Theoretical reaction models, which treat breakup as an excitation of the projectile to a 
two-body continuum state, most naturally express their results as cross sections describing 
the center of mass (cm.) and relative motions of the dissociated system, using 2-body 
kinematics. It has therefore been common for the experimental data to also be transformed 
to the cm. frame, for ease of comparison, e.g. the theoretical calculations of p0| , |21"| and the 
experimental data of JT7]]. This process is ambiguous in the case of inclusive data. Much 
more important is that the three-body cross sections are explicitly coherent in contributions 
from different spin-parity excitations of the projectile and so have the potential to offer a 
far greater insight into the projectile structure and the reaction mechanism. An excellent 
example of this is the interference observed |14| in the cross section of the 7 Be fragments, 
as a function of their component of momentum parallel to the beam direction, following 8 B 
breakup on a heavy target at 44 MeV/nucleon. 

In this paper we present calculations which are performed using full three-body kine- 
matics. These calculations are carried out within the framework of the coupled discretised 
continuum channels (CDCC) methodology, e.g. |]22"1J23|| , for breakup reactions of two-body 
projectiles. The interference between different excitation channels is shown to be signifi- 
cant for assessing the convergence of the calculations and those breakup excitations which 
contribute. The methods presented are applied to the breakup of 8 B on a 58 Ni target at 
Eiab — 25.8 MeV, for which new measurements have been reported |T^j24]1 . We compare the 
results of the full CDCC analysis, and also distorted waves Born approximation (DWBA) 
and truncated coupled channels calculations, with these available data for the laboratory 
angle and energy distributions of the 7 Be fragments. The calculations of Refs. ^I|^TJ showed 
the importance of higher-order breakup couplings, the couplings between continuum states, 
upon the 8 B* center of mass cross section angular distribution. We will show in this work that 
these higher-order effects are manifest even more significantly in the energy distributions of 
the 7 Be fragments following breakup. 
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II. THEORETICAL CONSIDERATIONS 



We consider the breakup reaction p — > c + v in which the projectile nucleus p is a bound 
state of a core particle c, of spin I and projection /i, and a valence particle of spin s 
and projection a. These particles are, presently, assumed structureless and so their internal 
wave functions are represented by the spinors Xj and X s . The total angular momentum 
of the ground state of p is J p , with projection M, the relative orbital angular momentum 
of the two constituents is £o, and their separation energy is So (> 0). The incident wave 
number of the projectile in the center of mass (cm.) frame of the projectile and target 
is K and the co-ordinate z-axis is chosen in the incident beam direction. The target t is 
assumed to have spin zero and no explicit target excitation is included. Target excitation is 
therefore present only through the complex effective interactions of c and v with the target. 
Our three-body solution of the Schrodinger equation calculates an approximate description 
of the projection of the full many-body p + t wave function onto the ground states of the 
target, core and valence nuclei. This three-body wave function is denoted ^fg M (r, R) where 

R is the position of the cm. of p relative to the target and f is the position of v relative to 
the core c. The particle masses are m p = m c + m v and m t . 



A. Construction of continuum bin states 



In the coupled discretised continuum channels (CDCC) method |g2|,^3l the breakup of 
p is assumed to populate a finite set of selected c + v excited configurations, with quantum 
numbers J p , £, j, where j = £ + s and J' p = j + I. Here, each of these spin-parity excitations 
will be assumed diagonal in all of these angular momentum labels. The excitations are also 
assumed to extend to some maximum relative energy S max (J') or wave number k max . This 
momentum range is then divided into a number N(J' p ) of intervals or bins, each with a width 
Aki = [hi — fci-i]. We label each such momentum bin by a = (i, J' p , £,j, s, I). 

In each of these relative motion bins a single representative wave function is constructed 
from those c + v scattering states f a (k, r) internal to the bin, with assumed angular momen- 
tum coupling 



YAr) ® X„ 



Xi 



J' V M< 



u a (r)/r . 



(1) 



The radial functions u a are square integrable and are a superposition 

/ 2 rki 



u a (r) 



7lN a Jki-i 



9a(k)f a (k,r)dk 



(2) 



of the scattering states, eigenstates of the c+v internal Hamiltonian H p , with weight function 



9a(k). 
for r - 



N a = 
-> oo, 



IkL \9a(k)\ 2 dk is a normalisation constant. The f a are defined here such that, 



fa(k, 



[cos 5 a (k)F e (kr) + sin S a {k)G i{kr)] 



(3) 



3 



where k G a and Fi and Ge are the regular and irregular partial wave Coulomb functions. So 
the f a are real when using a real c + v two-body interaction. An optimal discretisation of the 
continuum requires a consideration of the number, the boundaries k iy the widths Aki and 
the weights g a in the bins, which may depend on the J' configuration. Energy conservation 

relates the c + v cm. wave numbers K a and corresponding bin state energies £ a , as 

U 2 K a - Ti 2 Kl 

where we define each bin energy by £ a = (4> a \H p \4> a ) and where fx pt is the projectile-target 
reduced mass. 

For non-s-wave bins typically one uses g a (k) = 1 for a non-resonant continuum in which 
case Ni = Aki and = h 2 k 2 /(2/j, cv ) with kf = [kf — kf_ l ]/{3Aki). For s-wave bins it is 
an advantage to use g a (k) = k. This stabilises the extraction of the three-body transition 
amplitude at low relative breakup energies, discussed later in Eq. fli2|). In this case iVj = 
k 2 Aki and the bin energies are £{ = h 2 [kf — fcf_ 1 ]/(10/i ct; AA; i A; J 2 ). 



B. Coupled channels amplitudes 

These bin states <p a provide an orthonormal relative motion basis for the coupled channels 
solution of the three-body c + v + t wave function. The bins and the coupling potentials 
((p a \U(f, R)\4>p) are constructed, and the coupled equations are solved, using the coupled 
channels code FRESCO [E5|]. Here U(f, R) is the sum of the interactions of c and v with the 
target, which is expanded to a maximum specified multipole order q. The coupled equations 
solution generates the (two-body) scattering amplitudes, summed over partial waves, for 
populating each bin state J' , M' from initial state J p , M, as a function of the angle 6k of 
the cm. of the excited projectile in the cm. frame 



F M >M{K a ) = — J-Z (L0J p M\JM) (L'M - M'J'M'\JM) 
K V LVJ 

x exp(i[a L + a v \) \. S J LJp .. VJ ,{K a ) F L °(?o) Y L ^ M '(l a ) . (5) 

Here and oy are the Coulomb phases appropriate to the initial and final state cm. 
energies and the SLj p -.L'j^(K a ) are the partial wave S'-matrices for exciting bin state a with 
cm. wave number K a . When calculated using FRESCO |E3J, these amplitudes are expressed 



in a coordinate system with x-axis in the plane of K and K a , such that the azimuthal angle 
<pK a of K a is zero. When discussing three-body observables, it is more convenient to define 
the coordinate system with respect to the fixed positions of the detectors in the laboratory. 
For such a general x-coordinate axis the coupled channels amplitudes must subsequently be 
multiplied by exp(z[M — M']<Pk), with <pK referred to the chosen x-axis. 

For use in the following, the two-body inelastic amplitudes of Eq. (|5|) are re-normalised 
to that of the T-matrix by removal of their two-body phase space factors, so that 
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T M , M {K 



2-k% 2 Kr) 



fJ-pt V Ka 



^M'M{K a ) ■ 



(6) 



Throughout, we adopt scattering state and T-matrix normalisations such that, asymptoti- 
cally, the plane wave states exp(ik ■ r) that enter are multiplied by unity. 

It follows that the inelastic differential cross section angular distribution, in the center 
of mass frame, for excitation of a given bin state is 



da (a) 
dVL 



K 



2nh 2 



2J p + l 



2 ^e* I -a 



MM' 



(7) 



C. Three-body breakup amplitudes 

Less obvious is the relationship of the CDCC two-body inelastic amplitudes T M , M (K a ) 
to the breakup transition amplitudes T^-Mik, K) from an initial state J P ,M to a general 



physical three-body final state of the constituents • This is needed to make predictions 

for experiments with general detection geometries, since each detector configuration and 
detected fragment energy involves a distinct final state cm. wave vector K, breakup energy 
£ki and relative motion wave vector k. 

To clarify this connection we make the CDCC approximation to the exact (prior form) 
breakup transition amplitude, by replacing the exact c + v + 1 three-body wave function, 
M^f ■>")■> ^ its CDCC approximation ty CD , as 

K) = <0« (f) e iS -*\U(r, R)\^ M (r, R)} . (8) 

Here is the c + v final state. Upon inserting the set of all included bin-states, which 
are assumed complete within the model space used, then 

J?) = £ (4^|0f }(0f e^\U^^f oM {r,R)) , (9) 

at,M' 

where the sum is over all bins a which contain wave number k. We should now recognise 
that the matrix elements T M , M (K a ) of Eq. (Bp, obtained from the coupled channels solution, 
are precisely the transition matrix elements appearing in Eq. ([|), i.e. 

?M> M (K a ) = (0f e ig ^\U(r, mf oM {r, R)), (10) 

but calculated on the grid of 9 a and K a values. For the first term in Eq. (|) one obtains 

/^-)i^f ) = S)^J2HY(^sa\jm)(jmI^\J p M') exp[i5 a (k)] 9a(k)Yf(k) , (11) 



k ^ a 1 k^/N a , 
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where S a (k) = S a (k) + a a {k) is the sum of the nuclear and Coulomb phase shifts for c + v 
scattering at relative wave number k. It follows that the three-body breakup T-matrix can 
be written 

- - (2n) 3/2 

T^ :M (k, K) = {-if (W|jm) (jmI(i\f p M') exp[i5 a (k)} 

x Y e v (i.)g a (k)T M , M (a, K) . (12) 

Here the Tm>m{&, K) will be interpolated from the matrices T^, M (K a ), available on the 
calculated K a and 9x a grid. Specifically, for each value of K, we evaluate 

T M , M (a, K) = exp(i[M - M'](j> K ) \f^, M (K)/^ 1 (13) 

where the value of the bracketed term on the right hand side is interpolated from the coupled 
channels solution. 

In practice this interpolation is carried out as a function of the deviation of K from 
the threshold center of mass wave number. For non-s-wave breakup, the amplitudes are 
constrained to vanish at the breakup threshold K t h r , i.e. 

f> 2 K 2 f> 2 K 2 

W° M (Kt hr ) = 0, %^L = ^Jk-S . (14) 

We note that in Eqs. (|12"D and (0) the functional dependence of the T-matrix on the 
angles of k, the phase shifts S a (k), and the azimuthal angle 4>k are all treated exactly. 
The grid of 9k values can also be very fine without computing cost. The most important 
requirement is therefore that the number of bin states used to describe each [0 — > k max ] 
Jp excitation must be sufficient to allow an accurate interpolation of the amplitudes in the 
value of AK — \K — K t h r \, or alternatively in k. 



D. Three-body observables 



The three-body amplitudes of Eq. fli~2|) are used to compute the triple differential cross 
sections for breakup. If the energy of the core particle is measured then 



d 3 a 



2irfi pt 



1 



dVt c dVt v dE c h 2 K (2Jp + l) 



(15) 



With our T-matrix normalisations, and non-relativistic kinematics, the necessary three-body 
phase space factor p(E c ,Q c ,Q v ), the density of states per unit core particle energy interval 
for detection at solid angles Q v and fl c , is p7[ 



p(E c , fl c , fl v ) 



m c m v hk c hk v 

(27Tft) 6 



m v + m t + m v (k c - K tot ) ■ k v /k 2 



(16) 
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Here hk c and Kh v are the core and valence particle momenta in the final state and hK tot the 
total momentum of the system, all evaluated in the frame, cm. or laboratory, of interest. 
The association with the appropriate T-matrix elements in Eq. ( |T5| ) is made through 

K = k c + k v — — K tot , k = —k v - —k c . (17) 

m p + m t m p m p 



As the data under discussion here are inclusive with respect to the valence particle (pro- 
ton) angles, the calculated triple differential cross sections must be integrated numerically 
over Q v . The presented observables are also integrated and averaged over the solid angles 
AQ C subtended by the core particle detectors, with the stated detector efficiency profiles 
e(Q c ) |l7j], i.e. 



dtt c dE c / AQ C 



L da < ■ (18) 



It is most convenient to choose the x — z plane to be that defined by the beam and the core 
particle detector. 



III. APPLICATION TO SUB-COULOMB BREAKUP 



The method detailed above is applied to the breakup of B on Ni at energy Ei a b = 25.8 
MeV, for which new data are available [I7U23] . A first experiment was performed in 1996 
at the Nuclear Structure Laboratory of the University of Notre Dame (ND) ||10|| , one moti- 
vation being to clarify the importance of the E2 contribution to the Coulomb dissociation 
process, an issue which is still not completely resolved |I2"| . In that first experiment, the 
measured 7 Be fragments were detected at only one laboratory angle (~ 40°), assumed to be 
free from the influence of strong interaction contributions. However, as a result of theoretical 
predictions |^8[^9j of a strong nuclear peak beyond 40°, and claims also of Coulomb-nuclear 
interference at around 40°, a more complete experiment was recently carried out using the 
now upgraded ND facility. Measurements were obtained of an angular distribution of the 
7 Be fragments |IjJ and also of their energy distributions |24Tj for the range of measured 
laboratory angles. Although the removed proton is not observed, since the heavy fragment 
energies are identified, the presented 7 Be fragment distributions are known to contain no 
contribution from proton transfer reactions to bound states of 59 Cu. There may neverthe- 
less be contributions from knockout or stripping processes in which the proton excites the 
target and is absorbed. Such contributions are not calculated in this work. Proton transfer 
reactions to near-threshold (unbound) states of 59 Cu, if present, could also contribute. We 
comment briefly below on the latter. 



A. The CDCC model space 



Model space parameters common to all the CDCC calculations are as follows. Partial 
waves up to L max = 1000 and radii R up to R COU p = 500 fm were used for the computation 
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of the projectile-target relative motion wave functions. The continuum bins were calculated 
using radii r < r^n — 60 fm. The 7 Be intrinsic spin was neglected, the core being assumed 
to behave as a spectator. Thus we set 1 = 0. The proton spin, s = 1/2, was included and 
hence J' = j. 

In the final calculations presented all J' states consistent with relative orbital angular 
momenta t < 3, i.e. J' p up to /V/2, were included. We show that the effects of the g- 
wave continuum are small. The bin state discretisation was carried out up to maximum 
relative energy 8 max = 10 MeV for each state. The number of bins in the si/2-continuum 
was 32. For each of the other J p , 16 bins were used. These had equally spaced ki from 
k = to k max . In the case of the DWBA calculations shown the model space is the same, 
however, the bin states are coupled to the ground state in first-order only. Calculations 
using potential multipoles q < 4 in constructing the coupling potentials will be shown but 
the final calculations require q < 3. 

For the 7 Be- 58 Ni system, the interaction of Moroz et al. [^TJ] was used, as in the earlier 
analysis of Ref. fl2"0" |. The proton- 7 Be binding potential was taken from Esbensen and Bertsch 
(EB) . The model of Kim et al. |J2] is also considered. The potential used to construct 



the bin states was the same (real) potential as was used to bind the 8 B ground state, assumed 
a pure p%/2 proton single-particle state. The proton- 58 Ni potential is first taken from the 
global parameterization of Becchetti and Greenlees (BG) P3"f , but is also discussed below. 



B. Results of calculations 



It is important to note from the outset that the total breakup cross section angular 
distribution of the cm. of the excited projectile, the sum of the two-body inelastic differential 
cross sections of Eq. (0), is incoherent in the different bin components. This is not the case 
for the three-body amplitude of Eq. (|T2"D and the triple differential cross sections, Eq. (|15|). 
The practical convergence of the calculation, i.e. the dependence of the observables on the 
assumed model space, is therefore much more subtle in this case. 

The three-body calculations are found to require a more extended set of bins, excitation 
energies, and potential multipoles. Whereas the use of energy bins up to only 3 MeV 



of relative energy, and multipoles q < 2, e.g. in Ref. |20[| , gives stable (converged) cm. 
differential cross sections, in the sense of Eq. (0), this is not the case for the calculations of the 
triple differential cross sections and the energy and angle integrated distributions. We need 
the enlarged coupled channels model space, as detailed above, with bins extending beyond 
£max = 8 MeV to obtain a converged result for these three-body observables. Furthermore, 
even when the extended range of continuum energies is included, the bin discretisation may 
itself not be fine enough so that the basis of bin states is sufficiently complete. We have 
therefore verified the stability of our results, with regard to the bin size, by doubling the 
number of bins and confirming that the same results are produced. 
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1. Angular distributions 



The convergence of the three-body calculations with S max is clearly illustrated in Fig. [I]. 
Here we show the 7 Be laboratory differential cross section angular distributions from calcu- 
lations that include continuum bins up to £ max = 3,4,6,8, and 10 MeV. The calculations 
for this convergence test use multipoles q < 2 and £ < 3. The calculations use the BG 
proton-target potential and the EB proton- 7 Be potential. For the larger S max the bins have 
been constructed so as not to alter their low energy discretisation. The calculation of the 
three-body cross sections thus provides a different interpretation of the reaction mechanism, 
and evidence for significantly higher energy excitations than would be deduced from the 
earlier calculations and their comparison with the 8 B* cm. cross section. We will show 
that these high relative motion excitations are reflected in the calculated breakup energy 
distributions for 7 Be and the proton. 

Figs. |] and || present the calculated 7 Be laboratory differential cross section angular dis- 
tribution, integrated over energy and proton angles and averaged over the core detector solid 
angles, and compare this with the data [24]. The 7 Be detectors were circular, subtending 
a solid angle AQ C comprising a circle of radius 6° about the nominal laboratory angle 9i a b- 
They have a stated Gaussian efficiency profile e(9) with full width half maximum of 10.9° 
|T7| ]. Here 9 is measured from the nominal 9^ setting. 

The convergence of the calculations with multipole order, and also with the included 
continuum partial waves, is shown in Fig. 0. Here the long dashed curve is the result shown 
in Fig. [1], converged with respect to excitation energy, with q < 2 and £ < 3. The solid 
curve includes also the effects of the q = 3 multipole couplings for £ < 3. The dash-dot 
curve is a calculation where q = 4 multipole couplings and the £ = 4 breakup partial wave 
are included. The additional effects are small and the remaining calculations therefore use 
the truncated model space with q < 3 and £ < 3. 

The solid curve in Fig. |3] uses the BG proton-target potential and the EB proton- 7 Be 
potential. In Ref. [2~S| it was shown that different 7 Be- 58 Ni potential models give essentially 
the same shape for the 8 B* cm. angular distribution, while the cross section normalisation 
depends on the size of the 8 B g.s. wave function. The long-dashed curve in Fig. |3] shows 
the results of using the proton- 7 Be interaction of Kim et al. |32| . Consistent with earlier 
work, the cross section is enhanced due to the larger predicted 8 B root mean squared (r.m.s.) 
radius in this model. 

The Becchetti-Greenlees [[33| proton- 58 Ni potential, used above and previously, has sur- 
face imaginary strength and geometry parameters W=12 MeV, rv=1.32 fm, and 0^=0.534 
fm when computed at 3 MeV proton energy. Experience tells us [0 that the BG pa- 
rameters give reasonable fits to data only down to approximately 10 MeV. An alternative 
global parameterization, tailored for use below 20 MeV, has a similar imaginary strength but 
somewhat smaller radius and diffuseness parameters rv=1.25 fm, and avy=0.47 fm [55 and 
leads to very similar results. There are, however, also potential parameters fitted to elastic 
scattering data at 5.45 MeV J36],|34||. This analysis uses a Gaussian surface term and obtains 
a much reduced absorptive strength, 1^=3.5 MeV, r w =1.23 fm, and ay/=\.1 fm. We will 
refer to this as the VG potential. There is therefore some uncertainty in this potential input. 
The dot-dashed curve in Fig. |3| shows the calculated 7 Be angular distribution from the VG 
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potential. The cross section is changed only slightly at smaller angles. At the larger angles 
the calculated cross section is enhanced and is consistent with the experimental angular 
distribution data. 

Our calculations show that the 8 B structure (size) and proton-target potential uncer- 
tainties affect the calculations in characteristically different ways. The former produces an 
overall scaling while the latter produces, principally, a large angle enhancement. The data, 
currently, do not allow these effects to be discriminated further. In the final event, the 
overall agreement between the calculations and the data in Fig. |] is qualitatively similar to 
the comparisons made in Ref. [ 17| . There the calculated 8 B* cm. cross sections [^Oj .pl 
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compared with an approximate transformation of the measured 7 Be data of Fig. |3| to the 
cm. frame. Such approximate comparisons, however, are not necessary. 

We observe that the results of our calculations are qualitatively quite different to those 
presented in Ref. f37fl , where an isotropic approximation was assumed in calculating the 7 Be 
fragment laboratory cross sections. Those calculations show a radical change of shape of the 
angular distribution at forward angles which is not present in the calculations of Figs. [I], Inl- 
and [3] in which the angular dependences are treated exactly. 



2. Energy distributions 

In Fig. f| we show the calculated breakup energy distributions of the 7 Be fragments, 
together with the data from Ref. ||24|| , for four measured laboratory configurations. For the 
smallest angle, ~ 20°, the calculations and the data are the average of the distributions at 
Qiab = 19° and Q\ a b = 21°. For the largest angles, 50/60°, the curves and data are similarly 
the average of the distributions obtained at 9i a b = 50° and Q\ a b = 60°. The measured cross 
sections are zero outside of the range of the data points shown. The solid curves use the BG 
proton distortion and the EB proton- 7 Be potential. The general features of the data, their 
magnitude, centroids, and widths, are well described by the calculations. The long-dashed 
curves are the results using the Kim proton- 7 Be potential. They show an enhanced cross 
section, discussed earlier, but a very similar shape. The dot-dashed curves are calculated 
using the VG proton distortion and the EB proton- 7 Be potential. The small arrows on the 
energy axis in Fig. |] (and Fig. |^) indicate 7/8 of the 8 B energy for elastic scattering at each 
laboratory angle. An overall reduction in the mean energy of the heavy fragments within 
the breakup reaction is evident. 

Further insight is gained by looking at the results of DWBA calculations, and also cal- 
culations in which a subset of the continuum couplings are switched off, shown in Figs. 
||](a-d). The long-dashed lines show the DWBA calculations. The dot-dashed lines are the 
results of CDCC coupled channels calculations but in which all continuum-continuum (CC) 
couplings between bin states are removed. The solid lines are the full calculations, as were 
shown in Fig. |j. We see that the calculations in the absence of CC couplings, both DWBA 
and truncated coupled channels, show energy distributions which are strongly asymmetric 
and have an enhanced high energy peak. This is very similar to what is observed in the 
7 Be fragment parallel momentum distributions from 8 B breakup observed at higher-energy 



TM. As in that case, we show in Fig. that this asymmetry has its origin in the interfer- 
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ence between the El transitions to even breakup partial waves, and the E2 transitions to 
odd breakup partial waves. These EX amplitudes, which individually give approximately 
symmetric energy distributions, interfere to give strongly asymmetric responses. The very 
nearby kinematic cutoff in our case distorts the symmetry somewhat. The E2/E1 amplitude 
ratio in this lower energy case is also greater and so the asymmetry is enhanced compared 
to higher energies. 

In the full CDCC calculations these asymmetries are essentially removed as a result of 
the higher-order couplings. This higher-order coupling induced suppression of the E1/E2 
interference asymmetry was also a feature of the (higher energy) dynamical calculations in 



Ref. [31]. The suppression is more complete at the lower energy discussed here. Figure 
[7| shows the analogue of Fig. |(a), the calculated cross sections to odd and even breakup 
partial waves, from the full CDCC calculations using EB and BG potentials. Evident is 
the interference, both within and between the odd and even partial wave excitations. We 
note that the analogue of the E2 cross section, the p + / wave contribution, is not itself 
suppressed, and is in fact large. The interference between the two contributions in Figs. [7] 
and ^|(a) is however very different in the two cases. 

Also evident in these two figures is that the odd breakup partial waves contribution in the 
CDCC calculation is significantly narrower than that calculated using DWBA. This narrow- 
ing is already manifest in s +p wave two-step (q < 2 Coulomb) calculations and arises there 
from interference between the first-order E2 and second-order El amplitudes for populating 
the p wave continuum. The importance of these particular interfering paths was also noted 
in Ref. [31 , there in connection with a reduction in the calculated 8 B decay energy spectrum 



at higher energy, when going beyond first-order Coulomb excitation theory. The calculated 
energy distributions reveal even more clearly than those for the angular distribution the 
importance of a full treatment of the dynamical couplings within the continuum. 



3. Additional calculations and comments 



Since the proton separation energy from the 59 Cu(g.s.) is S p = 3.42 MeV, proton transfer 
to the 59 Cu(g.s.) would produce 7 Be fragments with « 26 MeV of kinetic energy in the cm. 
frame, and so such events are not part of the energy distributions measured. Those transfers 
that might contribute to the energy spectra of Fig. ^ would therefore be to excited (resonant) 
proton levels in 59 Cu* at around 9 MeV of excitation. If the proton- 58 Ni interaction sup- 
ported one or more potential resonances, then the CDCC reaction mechanism would include 
their dynamical effects since breakup, by projectile excitation and by transfer to unbound 
states, are not distinguishable mechanisms in the three-body reaction model used. Clearly, 
however, the ability of the proton- 58 Ni interaction to support such resonance strength, and 
its absorptive content, are closely related questions. As was noted earlier, Fig. [3|, use of the 
VG proton-target potential calculates an enhanced large angle cross section. Clarifying this 
sensitivity, and the possible role of such final-state resonances, requires further study and 
fine tuning of the proton-target potential. A full discussion of this topic is beyond the scope 
and motivation of the present article. 

With this sensitivity to the proton-target potential in mind, however, in Fig. ^ we 
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show the calculated proton laboratory angular distributions from the EB and Kim 8 B wave 
functions, and the BG and VG proton distorting potentials. We note that the magnitude, but 
not the shape, of the proton cross section angular distribution shows a significant sensitivity 
to the assumed absorption in the proton-target system. Precise data could therefore verify 
and constrain this element of the calculations. 

The shape of the calculated proton energy distribution, like that for the 7 Be fragments, 
shows little sensitivity to the absorptive content of the proton distortion or to the choice 
of 8 B binding potential. The calculations in Fig. |9| use the EB (solid) and Kim (long 
dashed) models for the proton- 7 Be interaction and the BG proton-target interaction. The 
dot-dashed curve uses the EB proton- 7 Be interaction and the VG proton-target interaction. 
The calculated proton energy distributions, integrated over all 7 Be fragment angles, peak 
for E p m 3.8 MeV and have a width r ~ 4 MeV. The tail of the energy distribution is seen 
to extend to high energy, reflecting the high relative energy excitations of the 8 B* discussed 
earlier in connection with the convergence of the CDCC calculations. Figure [10] shows the 
energy distributions predicted when the 7 Be fragments emerge at laboratory angles of 20, 
30 and 40°. In this case the arrows on the different curves indicate 1/8 of the 8 B energy 
for elastic scattering at each laboratory angle. The calculations show an increased average 
energy (acceleration) of the removed protons from the dynamics of the breakup process. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we have calculated the most exclusive three-body breakup observables of 
a two-body projectile using the coupled channels CDCC methodology. The formalism is 
applied to investigate the angular and energy distributions of the 7 Be fragments resulting 
from the sub-Coulomb breakup of 8 B on a 58 Ni target, the subject of recent experiments. 
We show that the convergence of the CDCC calculations of these observables is more subtle 
than that for the cross section of the cm. motion of the 8 B* and requires a significantly more 
extended space of 8 B* excitation energies. The required excitation energy range is clarified. 

Our calculations show that the 8 B structure and the absorptive content of the proton- 
target potentials affect the calculated 7 Be fragment angular distributions differently, the 
former producing an overall scaling, and the latter a large angle enhancement. Reducing 
the strength of the imaginary part of the proton potential in line with a phenomenological 
study |5fJ, provides agreement with the larger angle data. The full CDCC calculations are 
shown to provide a good description of the measured 7 Be fragment energy distributions. 
The widths and positions of these distributions are found to be rather insensitive to the 
details of the potentials used within the calculations. The presence of coupling between 
the continuum states is shown to be crucial to understand both the magnitudes of these 
energy distributions and their measured energy centroids. The absorptive content of the 
proton-target potentials affect the magnitudes of the calculated proton angular and energy 
distributions significantly, although their shapes are little affected. The calculated proton 
( 7 Be) fragment energy distribution reveals an overall increased (reduced) average energy of 
the fragment from the dynamics of the breakup process. 

The application of these techniques to calculate the parallel momentum distribution of 
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the heavy breakup fragments following the nuclear dissociation of the two-body system 11 Be 
will be reported elsewhere [38]. Further applications to systems with significant Coulomb 
dissociation strength, such as for 8 B breakup at energies of 40 MeV/nucleon and greater, 
are also in progress. 
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FIG. 1. Convergence of the calculated laboratory frame 7 Be cross section angular distribution 
following the breakup of 8 B on 58 Ni at 25.8 MeV as a function of the maximum proton- 7 Be relative 
energy included in the calculation. 
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FIG. 2. The calculated laboratory frame 7 Be cross section angular distribution following the 
breakup of 8 B on 58 Ni at 25.8 MeV. The long-dashed curve is the £ max = 10 MeV, I < 3, q < 2, 
calculation from Fig. [l]. The solid curve includes q = 3 multipole terms while the dot-dashed curve 
includes both q = 4 and I = 4 effects. 
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FIG. 3. The calculated laboratory frame 7 Be cross section angular distribution following the 
breakup of 8 B on 58 Ni at 25.8 MeV from the EB (solid) and Kim (dashed) models for the proton- 7 Be 
interaction and the BG proton-target interaction. The dot-dashed curve uses the EB proton- 7 Be 



interaction and the VG proton-target interaction. The experimental data are from Ref. [17 



17 




16 18 20 22 
E la /Be) (MeV) 



16 18 20 22 
E lab ( 7 Be) (MeV) 



FIG. 4. Calculated laboratory frame 7 Be cross section energy distributions following the 
breakup of 8 B on 58 Ni at 25.8 MeV for the laboratory angles indicated. The calculations use 
the EB (solid) and Kim (dashed) models for the proton- 7 Be interaction and the BG proton-target 
interaction. The dot-dashed curves use the EB proton- 7 Be interaction and the VG proton-target 
interaction. The arrows on the energy axis indicate 7/8 of the 8 B energy for elastic scattering at 
each laboratory angle. The experimental data are from Ref. [24|. 
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FIG. 5. Calculated laboratory frame Be cross section energy distributions following the 
breakup of 8 B on 58 Ni at 25.8 MeV for the laboratory angles indicated. The curves compare 
the full CDCC (solid), the CDCC in the absence of the CC bin couplings (dot-dashed), and the 
DWBA (long-dashed) calculations. All calculations use the EB 8 B ground state structure model 
and the BG proton distortion. The arrows on the energy axis indicate 7/8 of the 8 B energy for 
elastic scattering at each laboratory angle. 
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FIG. 6. Calculated laboratory frame 7 Be cross section energy distributions following the 
breakup of 8 B on 58 Ni at 25.8 MeV for the laboratory angles indicated. The curves show the 
separate odd and even breakup partial waves cross sections and their interference within the full 
DWBA calculation. 
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FIG. 7. Calculated laboratory frame 7 Be cross section energy distributions following the 
breakup of 8 B on 58 Ni at 25.8 MeV for the laboratory angle indicated. The curves show the 
odd and even breakup partial waves cross sections and their interference within the full CDCC 
calculation. 
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FIG. 8. Calculated laboratory frame proton cross section angular distributions following the 
breakup of 8 B on 58 Ni at 25.8 MeV, showing the role of the interaction between the proton and 
the target. The calculations use the EB (solid) and Kim (long dashed) models for the proton- 7 Be 
interaction and the BG proton-target interaction. The dot-dashed curve uses the EB proton- 7 Be 
interaction and the VG proton-target interaction. 
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FIG. 9. Calculated laboratory frame angle integrated proton cross section energy distributions 
following the breakup of 8 B on 58 Ni at 25.8 MeV. The calculations use the EB (solid) and Kim 
(long dashed) models for the proton- 7 Be interaction and the BG proton-target interaction. The 
dot-dashed curve uses the EB proton- 7 Be interaction and the VG proton-target interaction. 
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FIG. 10. Calculated laboratory frame proton cross section energy distributions following the 
breakup of 8 B on 58 Ni at 25.8 MeV for the 7 Be fragment laboratory angles indicated. The cal- 
culations use the EB proton- 7 Be interaction and the BG proton-target interaction. The arrows 
indicate 1/8 of the 8 B energy for elastic scattering at each laboratory angle. 
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